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J. Harer and D. Zagier have found a strikingly simple generating function [1, 2] for exact (all-genera) 1-point correlators in 
the Gaussian Hermitian matrix model. In this paper we generalize their result to 2-point correlators, using Toda integrability of 
the model. Remarkably, this exact 2-point correlation function turns out to be an elementary function - arctangent. Relation to 
the standard 2-point resolvents is pointed out. Some attempts of generalization to 3-point and higher functions are described. 

1 Introduction 

In quantum field theory, exact computation of correlation functions in all orders of perturbation theory is 
rarely possible. At best, we are able to find a few first terms, and study their properties. Only in low- 
dimensional and/or topological models, exact correlation functions can be sometime calculated. In this 
paper, we do such calculation in the Gaussian Hermitian matrix model [3]- [7], where the m-point correlators 
are given by the Gaussian integrals 
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over the space of N x N Hermitian matrices with flat measure, normalised so that (1) = 1. Originally 
designed to study random matrices, Hermitian matrix model is deeply connected to random surfaces [8]. 
This is because correlators Ci i ...i m (A r ) are polynomials in N with integer coefficients, like 
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C, A-. C* 4 = 2A 3 + A, C 6 r,.Y ; • 10.V-. C A A. C 2 .i = N 4 + 2N 2 , ... 

which count the number of matrix Feynman-t'Hooft graphs (ribbon graphs, fat graphs) of different genus, 
made of m vertices with i\, . . . ,i m legs. Such decomposition of correlators into several parts, related to sur- 
faces of different genera, is usually called genus decomposition. For example, correlator C4 gets contributions 
from two graphs with topology of a sphere (genus 0), and one graph of torus topology (genus 1): 
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Figure 1: Genus decomposition of the correlator (tr 4 ) = 2N 3 + N 



This connection to random surfaces makes matrix models applicable to a wide range of topics in modern 
physics, including two-dimensional quantum gravity [9, 10] and topological string theory [11]. Hcrmitian 
matrix model has a special place in this spectrum, being one of the simplest and most fundamental models 
(see [5] for its place in "M-thcory of matrix models"). 

It is therefore interesting and important to calculate Hcrmitian correlators. As usual in string theory 
[12], to represent the answers in a sensible form, it is more convenient to consider a generating function 
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E C il ... im {N)x\^...x% 

ii...i m =0 

However, these series do not converge, because the number of fat graphs of arbitrary genus grows too fast. 
If one wants to cure this problem, one can consider another series, with different normalisation of terms: 
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Historically, the first example of such a correlation function was constructed by Harer and Zagier [1] in the 
1-point case. They made the series convergent by special choice of double- factorial weight n(i): 



In fact, this choice of weight is quite natural: if one calculates the correlator C2i(N) = ^tr 2 *) with the 
help of Wick theorem, the total number of Wick pairings is exactly (2i — 1)!!. The result (2) becomes even 
simpler, if a generating function with respect to the matrix size N is also calculated [3]: 

00 x 2i \ N X 1 

E q ^ W (2iTT)fi = I-A (1 - A) - (1 + X)x* (3) 

It should be emphasized that, from the point of view of matrix model theory, this is a highly non-trivial 
relation: a generating function for 1-point correlators at all genera appears to be rational. In this paper 
we show, that similar relations can be established for the 2-point correlators: a clever choice of the weight 
n(ii,i2) makes the series convergent and, moreover, elementary functions. When both i\ and %2 are odd 
numbers, our result is the following: 
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' (27+ l)!!(2j + 1)!! ~ (A - If/ 2 ^ X - 1 + (X + l)(x 2 + y 2 ) 



(4) 



A generating function for correlators with even i\ and i 2 is similar, though a little more lengthy, see (54) 
below. As one can see, the 2-point generalization of the Harer-Zagier 1-point function is still an elemen- 
tary function. It is an open question, whether generalized Harer-Zagier correlation functions are always 
elementary. Affirmative answer to this question would greatly increase our understanding of the model. 
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It would be certainly interesting to generalize (3) and (4) to non-Gaussian matrix models, i.e, to Hermitian 
integrals with non-Gaussian weight. Especially interesting would be generalization to non-Gaussian models 
with multi-cut support [13, 14], which were recently related to supersymmetric gauge theories [15]. 

Another interesting direction of generalization are exact correlation functions in the presence of external 
field (matrix) 'J, see s.2 of [16]. We do not include this topic into the present paper, it will be discussed 
elsewhere [17]. The 1-point external-field correlation function turns out to be a simple deformation of the 
original answer (2): 



It is again interesting, that no special functions arise even in the external field case. We conclude, that 
non-standard correlation functions, especially of Harer-Zagier type (with double-factorial weights) appear to 
be simpler than the standard resolvents (with unity weights) and often can be calculated exactly. 

2 Correlators and integrability 

2.1 Generalities 

Before doing any actual calculations, let us briefly review the relevant properties of Hermitian matrix model. 
Partition function of the Hermitian matrix model depends on infinitely many variables tk known either as 
coupling constants or as time- variables. Partition function is a formal series in these variables 



where coefficients C^...^ are called m-point correlators. Instead of the full correlators (coefficients of Z^) 
one can consider connected correlators K il ,,, im , coefficients of the free energy Fn = logZjy: 



In terms of Feynman-t'Hooft diagrams (fat graphs) full correlators count all diagrams, while connected 
correlators count connected diagrams. This relation between the partition function and its logarithm is not 
specific for this model, it is a universal property of quantum field theory. At this point it is worth to mention 
a more traditional notation for correlators, coming from statistical physics: 




(5) 
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F N (t 1 ,t 2 , ...) = IogZjv(ti,t 2 , . . .) = Ki(N)ti + ^KijiN^itj + ^-.K ijk {N)titjt k + ... 




(6) 




(7) 



As a consequence of relation Fn = logZ^, the full and connected correlators are related by 
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Ki =Ci (8) 
Kij =C ij -C i C j (9) 
Kijk = Cijk CijCk CikCj CjkCi + ICiCjCk (10) 



and so on. Because of reflection symmetry of the action tr <j> 2 , correlators Ci 1 ...i m or Ki x _,,i m are non- vanishing 
only if i\ + . . . + i m is even. Note, that for N — 1 full correlators are quite simple: 

+?° ( - 1)!!, T,i = even 

Q,..i m (l)= / d4>^ e-^ l2 =\ (11) 
-i I 0, Sz- odd 

where £i = ii + . . . + i m is the sum of indices. To avoid confusion, we emphasize once again that in this 
paper we use normalised Gaussian integrals, i.e, the average of unity is unity: 

16 e-^l 2 = 1 
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Note also, that there is no difference between 1-point full and connected correlators: K = C%. The same is 
true for 2-point connected correlators with both odd indices: K 2 i+i,2j+i = C 2 i+i,2j+i, as a corollary of above 
identities and vanishing of 1-point correlators with odd indices, C 2 i+i = 0. Generally, connected correlators 
are somewhat simpler and we consider only them from now on. 

2.2 Virasoro constraints 

Usually evaluation of correlators in matrix models is done with the help of the loop equations, also known 
as Ward identities or Virasoro constraints [18]. In terms of the partition function, they can be written as 

dt b ^ dt a+b - 2 i+ ^_ 2 dtidtj 

In terms of free energy Fn — logZ^v, they take form 

9 v t 9Fn \- ST dp NdF N 

■jrr F N = >at a — + > — + > -7- — 7— , b>0 (13) 

a=0 i+j=b— 2 J i+]=b—2 J 

The technique based on loop equations [6, 18] allows to calculate the correlation functions for arbitrary 
genus. However, as also mentioned in ref.[3], it turns out that all-genera correlation functions, at least in 
the Gaussian phase, are much simpler deduced by another method, making explicit use of integrability of 
the model. A drawback of this method is that it is more difficult to generalize to non-Gaussian phases than 
the loop equation approach, but in the present paper we are only interested in Gaussian correlators. 
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2.3 Integrability 

Our main point in this paper is that correlators in the Hcrmitian matrix model are constrained by integrablc 
differential equations [7]. In principle, this should allow to calculate exactly all the quantities of interest. In 
the case of Gaussian Hermitian matrix model, the relevant equations are Toda equations [19], which can be 
written in terms of the partition function: 

Z N+1 Z N ^ 1 1 d 2 

Z% =NW 1 l ° gZN (W) 

Equivalently, these equations can be expressed in terms of free energy: 



F N+1 - 2F N + F N _! = log ( -t-^Fn ) (15) 



To calculate the correlation functions, we need to rewrite the Toda equations as relations between correlators. 
Differentiating eq. (15) by t i± . . . t im and using the connection between partition function and correlators 



one obtains the desired relations. Let us derive them explicitly, in the case of one- and two-point correlators. 
2.4 1-point correlators 

In the 1-point case, differentiating eq. (15) by ti we obtain 

Ki{N + 1) - 2K t (N) + K t (N - 1) = ^K^(N) (18) 

where we have used a simple identity K\^\ — Cn = N. Another identity we will need is the following: 

K i , hl (N)=i(i-l)K i _ 2 {N) (19) 

This identity is a corollary Virasoro constraints (13). Indeed, the Virasoro constraint for b = 1 implies 
Differentiating by t\ and using the Virasoro constraint once again, we get 
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F N = N+Y, P&rh fH £ + ]T a(a - 



(21) 



which justifies (19). Substituting (19) into (18), wc obtain 



Ki(N + 1) - 2K t (N) + Ki(N - 1) = ^^^(iV) 



(22) 



This system of recursive relations uniquely determines all the one-point correlators K \ (N) by reducing them 
to the N = 1 correlators which are already quite trivial: 



K 2i (l) = C 2i (l) = y # e^ 2 - (2i - 1)!! 



(23) 



Interestingly, 1-point correlators satisfy additional recursive relations, different from (22): 



Ki(N + 1) - K t (N - 1) = l -±^Ki(N) 



(24) 



These relations are similar to (22), but look slightly simpler. However, they do not follow from Toda 
equations and their generalization to 2-point and higher-point correlators so far remains unclear. We do not 
use them in this paper, but it is important to mention that they exist. 

2.5 2-point correlators 

We now do the analogous calculation in the two-point case. Differentiating eq. (15) by U and tj, we obtain 
Ka(N + 1) - 2K tJ (N) + K^N - 1) = ±K iijAA (N) - ±K hlA (N)K jA AN) (25) 



As a corollary of Virasoro constraints, 

Kij,i,x(N) = i(i - l)Ki_ 2tj (N) + 2ijK i _ 1J _ 1 (N)+j(j - l)K itj _ 2 (N) 

Substituting (26) into (25), we obtain recursive relations for the 2-point correlators: 





-2K ij (N) + K ij (N-l) = - 1 -^ 


ffi 1) K i . 2 (N)K j . 2 (N) + 


1 

+ N 


U{i-l)K i _ 2 , j (N) + 2ijK i _ l!j _ 1 


N)+j(j-l)K id _ 2 (N)^ 



(26) 



(27) 
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In complete analogue with the 1-point case, these recursive relations allow to express arbitrary correlators 
Kij(N) through the N = 1 correlators, which are simple according to (11). Note however, that there are 
two essentially different cases: when both indices are even 

K^jW = C 2i;2j (l) - C 2i (l)C 2j (l) = (2i + 2j - 1)!! - (2i - l)!!(2j - 1)!! (28) 

and when both indices are odd 

tf 2 »+i,2i+i(l) = C 2i+ i, 2j+1 (l) - C 2i+1 (l)C 2j+ i(l) = (2i + 2j + 1)!! (29) 

The second "totally odd" case is somewhat simpler, because in this case the quadratic contribution the 
product of two odd 1-point correlators - vanishes. This is actually one of the reasons, why the totally odd 
generating function (4) is simpler, than its totally even counterpart. 

3 Correlation functions 

There are many different ways to solve discrete relations like (22) or (27). For example, one can start solving 
them iteratively and try to guess a combinatorial formula. As explained in [3] , a more elegant way to proceed 
is to pass from particular correlators (which depend on discrete indices like i or j) to generating functions 
(which depend on continuous variables). In terms of generating functions, the discrete relations like (22) 
turn into differential equations which often admit simple and elegant solutions. 

For the sake of brevity, we call generating functions for correlators simply correlation functions. These 
objects are also known in literature as (multi-)densities [3]. In this paper, we restrict consideration to three 
types of correlation functions. First of all, the standard correlation functions, with unity weights: 

PN{Xl ,..., Xm ) = ((tr ( T ^)...tr ( I ^)))= i ±j i ,.UN) X r... X t 00) 

These functions are the most commonly used in literature, especially in the context of loop equations 
approach [6]. To avoid misunderstanding, let us note that in the context of loop equations the same functions 
are usually defined in a slightly different way 

W N ( Xl ,..., Xm ) = ((tr (^)-tr (^)))= ii £ o ^-^W^- 1 ...^- 1 (3D 

and known as resolvents. We use the terms "standard correlation function" and "resolvent" as synonims 
here, because it is not a problem to transform one into another: two different definitions are related by 

Wn{xi, ■ ■ .,x m ) = pn ( — , — J 

X\ . . . X m \ X 1 x m J 

Therefore, it does not make a big difference. In this paper, we prefer to use functions pjq rather than 
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functions Wn- Second, we consider the exponential correlation functions, i.e, with factorial weight: 



00 11 i 

e N (x u ...,x m ) = ((tv (e^)...tr (e*™*) )) = V K^.^N) X }"- X ™ (32) 
\\ 11. %\\ . . . i m . 

Jl...J m =U 

The third functions are those which play the central role in our paper - the Harer-Zagier correlation functions: 

00 hi f (i — 1)!!, i — even 

<p N (xi,...,x m )= J2 K n-im(N) *\-- X f. v n(t)={ (33) 
i!...i m =o n^!j...n^ ro j ^ ., !; z = odd 

Obvoiusly, different choices of weights can be useful under different circumstances. The Harer-Zagier func- 
tions arc useful just because they provide simple and explicit answers: it is enough to take a look at (3) or 
(4). The exponential and standard correlation functions are useful for another reason: it is because Toda 
equations, rewritten as differential equations on these functions, take the simplest form. This can be seen 
already for the 1-point equations (22), since the term i(i — l)Ki- 2 suggests two natural choices of weights: 

and 

The first choice corresponds to the exponential correlation function, while the second choice corresponds to 
the resolvent in its usual form (31). For other choices of weights, operator in the right hand side would be 
more complicated than just x squared or the second derivative. That is why, in a scncc, these two choices 
are singled out by the equation itself. 

Of course, the three correlation functions are related by various integral transformations. First, the 
standard and exponential functions are related just by Laplace transform: 

00 00 

p N {x 1 ,...,x m ) = J dy 1 . . . J dy m e N {x 1 y u . . . ,x m y m ) e ~ Vl ~-~ Vm (34) 


Second, the standard and Harer-Zagier functions are related by a certain Gaussian transform: 



oc 00 

p N (xi, . • • , x m ) = J (1 + yi)dyi ... J (1 + y m )dy m <Pn{xiVi, x m y m ) e -f?/ 2 -----f™/ 2 (35) 



Third, the exponential and Harer-Zagier functions are related by a contour-integral transform: 



ipNiyi, ■ ■ ■ ,Um) ( x\ x 2 



e N (xi, x m ) = f (1 + yi)dyi ... f(l + y m )dy m '-^ exp I + . . . + (36) 

Vl ■ ■ ■ Dm \ ^Vm J 
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If one is able to find one of these functions - e, p or (p - it is not a problem to convert it into another. Therefore, 
one can freely change the weights in order to simplify the solution. Another important simplification, which 
was suggested in [3] and which we intensively use, is to consider the universal generating functions, i.e, 
generating functions w.r.t N with parameter A: 



p(X;x 1 ,...,x m ) = ^ N Pn(xi, ■ ■ ■ ,x m ) (37) 

N=0 



e{X\x 1 ,...,x m ) = ^ X N e N (xi, . . . , x m ) (38) 

N=0 



ip(X;xi, . . .,x m ) = y X N (p N (xi,. . . ,x m ) (39) 

As we will see below, transition to universal functions greatly simplifies both the equations and the answer. 
We are now going to rewrite eqs. (22) and (27) as differential equations on correlation functions and solve 
them directly. The solution in the 1-point case is the well-known Harer-Zagier function (3). In the 2-point 
case, its generalization is obtained. 

4 Harer-Zagier correlation functions 
4.1 1-point function 

Rewritten in terms of the generating function for 1-point correlators 



cq. (22) becomes a differential equation of first order: 

1 d 



ifiN+i(x) - 2ip N (x) + Lp N -i{x) = —x—(x 2 Lp N {x)^ (40) 



Passing to generating functions with respect to N, we obtain 

( ^^(A; x) \ = x^- (zV(A; x )) (41) 



dX V A J dx 



This is a first-order partial differential equation with general solution 



The function U (z) is determined from the initital condition 
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|^(A;*) 



+ 00 

00 2i r 00 1 

o = ^ =l( «) = £ / # ^ e -^ 2 / 2 = = _J 

i=0 v _ i=0 



(43) 



In other words, the initial condition is merely a consequence of the fact that Ki(l) = (2i— 1)!! and the choice 
of weights. Comparing the general solution with the initial condition, we obtain U(z) = 1/(1 + z) and 



(44) 



In this way the Harer-Zagier correlation function </?(A; x) can be found as solution to a linear differential 
equation of first order in two variables A and x. Note, that the universal A-dependent correlation function, 
which contains information for all dimensions N, seems to be very similar to the N = 1 function: namely, 
they are related by multiplicative transform 




2 1 + A 2 A 

x ^— X > ^(i~Ay 



(45) 



Interestingly, this property has a literal analogue in the 2-point case, see (55) below. 
4.2 2-point function 

For two-point correlators, the two separate generating functions can be introduced, even and odd: 



h3=0 



(46) 



V~ N {x,y) = E^W^( JV ) (2 i + i ) !!( 2 j + l)!! 



h3=0 



(47) 



In terms of these generating functions, eq. (27) becomes a system of two equations: 



<pn+i( x > v) - 2i pn( x ' v) + <pn-i( x > y) 



1 ( d 2 d 2 \ , . . 2 d 2 _. , 1 / d 2 2 2 \ , n / n 
N \ X Tx X + W )^ y)+ N Xy dx^ Xy ^ ix ' y) ~N^ [xy—x y ) v N (x)<p N (y) (48) 



and 



<P N+ i{x,y) - 2ip N (x,y) + (p N _ 1 (x,y) = (x 2 -^x + J/ 2 J^ t P N (x,y) + ^: xy Lp+{x,y) (49) 



Passing to generating functions with respect to N, we obtain 
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f x 9 A 1 -A) 2 +a 



d 



dx a dy' 



i) 2 



+ y~x~y 2 )^ + ( A ! ' T ' V) + "a - H -2 ^ (^ £, 2/) ~ G(A; ir, 2/) 



dxdy 



A— 



a /(i-A) 2 



<9A V A 
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v (A;z,t/) ) = ( x2 q^ x + y 2 Q^y ) v ( x ; x >y) + 2 %y<p + (*;x,y) 



(50) 



or in a matrix form 



/ A 9 

<9A A 



d 9 9 9 



n 52 

-2xy———xy 



dxdy 



d (l-A) 2 
SA A 



9 9 9 a 
x'—x-y'—y . 
ox ay ) 



G 




where the free term G(A; x, y) is given by 



G(A; :r, y) = ^ — [x y ^-^- x 2 y 2 j tp N (x) <p N (y) = A 



2xy 



(\-l)(l + x 2 y 2 ) + (\ + l){x 2 +y 2 ) 



(51) 



This system of two differential equations is more complicated, than in the 1-point case, but it is still linear 
and can be solved by elementary means. First of all, we need to find the initial conditions, what is done 
again by comparing to the N = 1 case. For N = 1, we know the correlators explicitly 

K 2i ,M = Cm,!«(1) - C 2i {l)C 2j {\) - (2i + 2j - 1)!! - (2* - l)!!(2j - 1)!! 
^+1,2^+1(1) = C 2i+1 , 2j - +1 (l) = (2t + 2j + 1)!! 
so we can compute the generating functions 



<Pn=i(*, y) - 2^ (2i + l)!!( 2 j + 1)!! X y 



: arctan 



^/l — x 2 — y 2 \ \J\ — x 2 — y 2 



v + N=1 {x,y) = 

i,3=0 



* + 23 - 1>!! - * - 1,!!,2J - .^2,(4-4, 



(2i-l)!!(2j-l)!! 



x 2 — y 2 \ dx dy 



Therefore, initial conditions for (50) are 



arctan 



^2/ 



x — 2/ 
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— ^ + (A;x,?/) =^=i(a;,2/)= ., . 
o\ a=o x A — y z \ ox Oy 



xy l d d , _ 

r AT ViV=l(*»tf) 



(52) 



The unique solution of (50) with initial conditions (52) is the following: 
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ip (X;x,y) 



( xy^/X - 1 

arctan — . 
A \^X-l + (X + l)(x 2 +y 2 ) / 

(A - l) 3 /2 v /A-l + (A + l)(x2 + y2) 



xy i Q d 

x a V7T ) <P~(^ x ,y) 



X 2 — J/ 2 y (fo Qy 



(53) 



(54) 



A(A + l)x 2 y 2 
(1-A) 



(A-l + (l + A)(.T 2 +y 2 )) 1 (A-l + (l + A)( a ; 2 + 2/ 2 ) + (A-l) a ;V) 



A(A + ljsg / 1 



-3/2 



xy^/X — 1 



— 4ttt A - 1 + (A + l)(x 2 + y 2 )) arctan — ; = 

(1-A) 3 M 1 A '"V I v /A-l + (A + l)(.x 2 +y 2 ) 



It is an elementary exercise to substitute (53) and (54) into (50) and check that equations are satisfied. The 
check of initial conditions (52) is equally simple. Just like in the 1-point case, the universal A-dependent 
correlation function is related to the N = 1 correlation function by a simple transformation 



2 9 1 A . o 2\ A 

- -{x + y), xy^xy, ip ^ — -tp 



* +y 



1-A 1 



(1-A) 2 



(55) 



This striking relation, valid for both even and odd functions, is clearly a hint for some larger structure, which 
can be completely revealed only by the study of the 3-point and higher-point cases. Note also, that it is 
straightforward to extract correlation functions for particular TV from the universal A-dependent correlation 
function. However, the formulas become less explicit, i.e, involving an integral: 



<Pn 



xy 

(x, y) = j 7 



dt 



l + t 2 + x 2 + y 2 ^ A 



2(x 2 +y 2 ) \\l + t 2 -x 2 -y< 



(56) 



+ 7 xydt ( ( l + t 2 + x 2 + y 2 \ N 2N(x 2 + y 2 )(l+t 2 ) + (x 2 +y 2 ) 2 -(l + t 2 ) 2 \ 

<PnK*>V) J ( x 2 +y 2)2 \y i + t 2 _ x 2 _ y 2) (l + t 2 + x 2 +y 2^ 1 + t 2_ x 2_ y 2} + J 



(57) 



Formulas (53) and (54) completely describe the exact 2-point correlators: it is enough to write 



tr ^ +1 tr (t> 2l+1 )) 
(2i + l)!!(2j+l)!! 



arctan 



xy^JX — 1 



coefficient of x 2k+1 y 2m+1 X N in 



v /A-l + (A + l)(x 2 +jy 2 ) 



(A-l) 3 / 2 ^X-l + (X + l)(x 2 +y 2 ) 



and similarly for even correlators. Generalisation of these formulas to 3-point and higher-point cases is not 
straightforward, since Toda equations become more complicated and explicit solution becomes increasingly 
hard to find. Generalization to non-Gaussian (say, Djkgraaf-Vafa [14]) models is even more obscure, since 
integrable equations (14) are non-trivially modified in these models. 
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5 Exponential correlation functions 

5.1 Recursive relations 

Direct generalization of (53) and (54) looks problematic. To bypass those difficulties, let us use the freedom 
in the choice of weight and consider not Harer-Zagier but exponential correlation functions: 

oo i* 

e N (x 1 ,...,x m ) = ((tv (e^)...tr (e*™*) )) = £ K n ... lm (N) ^ 



%\\ . . . i m ! 

«i ...t m — o 



Expressed in terms of these functions, the Toda equations do not contain any differential operators at all: 



x 2 



e N+ i(x) + e N -i(x) = 2e N (x) + — e N (x) 



frg _|_ ^,2 y2 

e N +i{x,y) + e N ^ 1 {x,y) = 2e N (x,y) + — e N (x,y) - ——e N (x)e N (y) 



(j, _|_ y _|_ (x ~\- y)^ z^ 
e N+1 (x, y, z) + e N -i(x, y, z) = 2e N (x 1 y, z) H — e N (x, y, z) — ~^e N {x, y)e N {z)- 

_|_ y^ (y -\- z) X^ y^ Z^" 

jf — j^e N (x, z)e N (y) — — jje N (y, z)e N (x) + 2— — —e N (x)e N (y)e N (z) 



and so on, generally 



eN+i + eN-i = c N e N + g N (58) 



where 

2 

c N = 2 + —{xi + . . . + x m f 

and gN(xi, ■ . ■ , x m ) is the function which can be considered as already known — by recursion: 

gN(x) = 

2 2 

x y 

9N(x,y) = - — — e N {x)e N (y) 

(x -h y) z^ (x -h Z} y^ 
g N (x,y,z) = — — e N (x,y)e N (z) — — — e N (x, z)e N (y)- 

(y _|_ x^ x^ y^ z^ 
jj — e ^(y, z)e N (x) + 2— — —e N (x)e N (y)e N (z) 

We now turn to explicit solution of equations (58). In contrast with the Harer-Zagier case, where solution of 
Toda equations is a non-trivial procedure, in the exponential case solution has nothing to do with differential 
equations and is rather simple. We emphasise, that this simplicity is due to the choice of weight. 
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5.2 Determinantal solution 

As usual for integrable equations, solution of cq. (58) can be given explicitly in terms of N x TV determinants: 



cm = det 

NxN 



( 


ei 


-5i 


52 


-53 • 


.. \ 






1 


Cl 


1 


. 











1 




1 




(59) 










1 


C3 






V 










•■ / 





To prove this, it suffices to expand (59) by elements of the last two rows. For example, for N — 3 we have 



e3 



ei -51 52 
1 ci 1 
1 c 2 



C2 



ei 


-5i 




ei 


52 


1 


Cl 




1 


1 



C2e 2 + 52 - ei 



so that (58) holds, for A = 4 we have 



e 4 



ei 


-5i 


1 


Cl 





1 









(•?, 



ei 


-0i 


52 




ei 


-0i 


-53 






1 


Cl 


1 




1 


Cl 





= c 3 e 3 + g 3 - e 2 


(61) 





1 


C2 







1 


1 







1 

C2 1 
1 C 3 



and (58) holds again. Generalisation is obvious. 
5.3 Orthogonal polynomials 

If, instead, one expands the determinant by elements of the first row, one obtains 

• \ 



ex = ei det 



/ ci 1 

1 C 2 1 
1 C 3 

v 



N-l 



+ ^2 9i det 



»=i 



/ Q+i 1 

1 Cj+2 1 
1 C i+3 

v 



... \ 



■ / 



Therefore, the answer for the arbitrary m-point exponential correlation function can be written as 



JV-l 

e N (xi, . . .,x m ) = T^(Si) 2 j ei(xi, . . .,x m ) + ^ T^^x) 2 ^ gi{x u . . .,x m ) 



(60) 



(62) 



(63) 
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where TUu) are special polynomials 



/ 2 



!?(«) = 



det 

0-i) x 



+ 



2+1 



2 + - 



i + 2 
1 2 



i + 3 



(64) 



with a general formula 



0-i)/2 a u , 

T >) = E E (-i)°n( 2 - 

a=0 fci+...+feL=0 1=1 ^ 



i + / + 2ki + . . . + 2ki 



L = j - 2a - 1 (65) 



A few first polynomials TUu) are: 



T » = 0, ^ = 0, T 2 = 0, ... 
!?(«) = 1, ^ = 1, T? = l, ... 

T 2 °=7i + 2, T 2 1 = i U + 2, T 2 2 = i,i + 2, ... 
T 3 ° = i u 2 + 3 U + 3, T^ = ^u 2 + ^u + 3, Ti = ^u 2 + 7 -u + 3, 



u - 6 « 



3 , 43 
' ' -24 U " + 4 tt : 12 



2u 2 + 6u + 4, T| = — u 3 + ^u 2 + ~~u + 4, If 



1 3 2 2 13 



It is easy to see, that polynomials TUu) satisfy recursive relations 



i + j 



,T j (u)=T j+1 (u)-2T](u)+T]_ 1 (u) 



(66) 



which can be viewed as three-term relations in orthogonal polynomial theory. The three-term relation implies, 
that polynomials T*(u) form an orthogonal set of polynomials with respect to some, yet unidentified, local 
measure [7]. For example, for i = they are orthogonal on the segment (— oo,0) with measure ue u du: 



u 

/ 



T?(u)T°(u)ue u du = j5 jk 



(67) 



i.e, for i = they belong to the family of generalized Laguerre polynomials: T^(u) = Laguerre}^ ( — u). It 
would be interesting to find the local measure, corresponding to polynomials Tj(u) for arbitrary i > (its 
existence is a consequence of three-term relations). In the next section, we derive an integral equation on it. 
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5.4 The local measure 

For convenience, let us introduce normalised (with unit leading coefficient) polynomials 

Q){u)= it+J ~ 1)l T;(u)=uJ- 1 + ... (68) 

which satisfy recursive relations 

uQ){u) = Q) +1 (u) 2{ l + ] )Q)(u) + ^±^Qj_ l(u ) (69) 
A first few polynomials Q l j{u) arc 

Qh(u) = 
Q\{u) = 1 
Q\{u) =u + 2i + 2 
Ql(u) = u 2 + (4i + 6)u + 3i 2 + 9i + 6 

and so on. For each particular value of i, the system of polynomials {Qj},j = 0, 1, 2, ... is orthogonal with 
respect to the unknown measure dpn{u) = uji(u)du. The orthogonality relations can be written as 

o 

J Q){u)Ql{u)u i {u)du = Q, j^k 

— oo 



(we assume that the segment is always (— oo,0), just like for i = 0). The moments Mfc(i) = J u k tOi(u)du 

— oo 

can be found from these orthogonality relations. Several first moments, obtained in this way, are 

M (i) = l, Mi 00 = -(2i + 2), M 2 (i) = 5i 2 + Hi + 6, M 3 (i) = -(Ui 3 + bOi 2 + 60i + 24), ... (70) 
Direct calculation shows, that generating function for these moments has a form 

oo 7 2 7 3 7 4 

E (-l) fc M fc — — ■ = z+(2i + 2)— + {bi 2 + Hi + 6)— + (14* 3 + bOi 2 + 60* + 24)— + . . . = 

k=0 K + i I & 4 

(71) 

1 . A .,. z(i + l) 2 (« + 2) 2 i(« + l) 2 (i + 2) 2 (i + 3) , \ 
= .,. , log 1 + i(t + l)z + -i ^ ^z 2 + -i y v y v >-z 3 + ...) 

This relation can be written as an integral equation 

?, , , 1 , + + l)\z k \ 1 n , , . 

7 log(1 + ^"V 2 * 4 = log QC iid-i)! ¥ J = ^TT) log 2Fo (*' 2 + 1; z) (72) 

where 2-F0 is the generalized hypergeometric function. The local measure Ui{u)du can be found as its solution. 
It in present paper, we do not solve this equation and we do not use this measure in what follows. 
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5.5 1,2,3-point functions for particular N 

Application of eq. (63) to particular correlation functions is straightforward. For m ~ 1 we have 



e N 



(x)=T° N (x 2 )e 1 (x) (73) 



for m = 2 we have 



e N (x,y) = T° f ((x + y) 2 )e 1 (x,y)- £ ((a: + y) 2 )^)^) (74) 



for m = 3 we have 



JV-1 



ejv (x,y,z) - T°((x + y + z) 2 ) ei (x,y,z) + 2^ TAr((a ; + y + z) 2 )^^e,(x) e4 (y) e4 (z)- (75) 



£ T ^ ((* + V + z ?) ( {XZ % yZ? e t (x + y)e l {z) + {xy ^ e t {x + z)e i {y) + {xy ^ e t {y + z)e i {x)\ 

i=l ^ ' 

and so on. In this way, all the exponential correlation functions ej^(x\, . . . ,x m ) are expressed through 
polynomials TJj(u). These expressions are fully explicit and constructive (since all kinds of explicit formulas 
are available for T ? (u)) but they are not very illuminating, certainly less beautiful than eq.(4). 

5.6 Universal 1,2,3-point functions 

To cure this problem, one may switch to universal A-dependent correlation functions, hoping that this will 
improve the situation. After conversion to universal functions 



e(A;zi,...,a; TO ) = £ X N e N (x 1} . . . , x m ) (76) 

AT=0 

the Toda equation (58) transforms into a differential equation of first order: 

A Ja~ ( ~ A~^ e(X;xi,...,x m ) - g(X;x 1 , . . . ,x m )j = (Sa;) 2 e(A; x ly . . . , x m ) (77) 
where T,x = x\ + . . . + x m and g(X; x\, . . . , x m ) is the differential equation's free term: 

oc 

g{X;x 1 ,...,x rn ) = xN 9n(xi, ■ ■ ■ ,x m ) (78) 
Unique solution of eq. (77), satisfying the usual N = 1 initial condition 
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^e(A; x\ , . . . , Xm) 



A=0 



— 61 (xi, . . . , x m *j 



(79) 



is given by 



e(A;a;) = (r3AF exp (T^rJ ei(a;) + (i-A)^ - ,-;/ 





(l-A)(l-t) 



(80) 



Application of this formula to particular correlation functions is straightforward. Let us begin with the 
simplest case, i.e, with the 1-point correlation function. Since g(X; x) = 0, it immediately follows, that 



e(X,x) = \ cxp ( X — I ei(x) 



(1 - xy 



1 - A 



Moreover, the initial function can be calculated explicitly: 



+ OC 



(81) 



so that 




(82) 



This answer is consistent with the Harer-Zagier 1-point function (3) and transformation (36): 



e(A; x) = <j> 



(i + y)dy 
y 



p(A; y) cxp ( ^ 



2 \ / 2 

x \ fx 



X r(l + y)dy eXP VV/ A rdy exp \2y 2 J X fx 2 1 + A 



(1-A) 2 / y }_A_ V 2 (1-A) 2 7 y 1-A (1-A) 2 ° XP U 1-aJ (83) 

l+A y 1+A y 

The last equality is a direct corollary of Cauchy residue theorem. Let us turn to the ncxt-to-simplest case, 
i.e, to the exponential 2-point function. In this case, the free term g(X; x, y) is no longer zero: 

00 x 2 2 
g(X;x,y) = - ^ X N -j^-e N {x)e N {y) = 

N=0 



du 1 du 2 ( v-^ —N —N\n\ ( x 2 1 + , y 2 l + u 2 \ 
— (J> « 2 A Jex P ^ y . T — - + y . T —-j 



duidu 2 fx 2 l + ii! y 2 l + u 2 \ 

' ex P 1 "TT ' 1 + "TT ' 1 ( 84 ) 



A — uiu 2 V 2 1 — m 2 l — u 2 
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Consequently, the 2-point function can be written as 



e(A; x, y) = exp ( ^ ^~ V \ X ) a (x, y) - 



dt 



A 



du\d,U2 



■ exp 



1 - A 



(x + y) 2 (\-t) x 2 l + m y 2 l + u 2 



{t- Ul u 2 ) 2 ± V( 1 - A )( 1 -<) 2 2 l-u 5 



(85) 



where 



ei {x,y) = J ^ e -^ 2 /2+^+y0_ j y Q 



d(j>e 



-<P 2 /2+x4 



d<fie 



2 /2+y4, I = e (x+s/) 2 /2 _ e x 2 /2 e y 2 /2 



In principle, it should be possible to obtain the same formula with two contour integrals in a different way, 
directly applying the transformation (36) to the Harer-Zagier 2-point function. Note, that (85) is far less 
concise, than (4). Instead, its advantage is the possibility of generalization: it is not a problem to write its 
analogue for any m-point function. Say, for the 3-point function we have 



A 

/ 



dt 



(1-A) 2 . f (x + y + z) 2 \ \ 
^ e{X;x,y,z) = exp I — J ei(x,y,z) 



(y + z) 2 du\du 2 ( (x + y + z) 2 (A - t) x 2 1 + m 



t(t - U1U2) 



■ exp 



(l-A)(l-i) + 2 1-m 



e(u 2 ;y,z) 



A 

+ /* 





(x + z) 2 du\du2 
t(t — U1U2) 



exp 



(x + y + z) 2 (X-t) y 2 l + ui \ 



+ 



dt 



A 



(x + y) 2 du 1 du 2 f(x + y + z) 2 (X-t) z 2 l + m\ 

ex P I — m — ~ 1" — ' — ) e \ u 2]x,y) + 



t(t - U1U2) 



(l-A)(l-i) 2 1-m 



2du\du 2 du 3 ( [x + y + z) 2 (\ - t) x 2 1 + u\ y 2 1 + U2 z 2 1 + u 3 
(t - u lU 2U 3 ) 2 6XP V (1 - ~t) + Tl - m + ~2~ 1 - U2 + ~2l-u 3 



(86) 



where 



ei(x, y, z) - e ( x+ » + *> / 2 - e (*4*)7V / 2 _ ^x+zY/^/i _ e (y+zY 12^/2 + 2e *W / 2 e* / 2 



Clearly, functions (82), (85) and (86) belong to a family of exact solutions, which are less elegant, than (4). 
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6 Standard correlation functions (resolvents) 
6.1 Genus expansion 

Among the all correlation functions, the most widely used ones are the standard correlation functions: 



PN ( Xl ,..., Xm ) = ((tr ( T J ^)...tr ( Y ^)))= i ±^-..UN)^...^ (87) 
As we already mentioned in s.3, they are better known as resolvents and usually defined as 

W N (x 1 ,...,x m ) = ((tr (^)-tr (^))>= t E^^W»r <1 - 1 -"^- 1 ( 88 ) 

All-genera resolvents are divergent series, because the number of fat graphs of arbitrary genus grows rapidly. 
This property is typical for any perturbation theory, where the number of Feynman diagrams of order n 
usually grows as nl. In practice this means that one can at best hope to represent the full pn as an integral, 
just like it happens with the archetypical divergent sum 

oo 

which is actually divergent only for x € R + . Such integral representation for resolvents is naturally provided 
by Harer-Zagier correlation functions: if the latter are known, the resolvents are given by 



oo oo 

p N (xi, . . .,x m ) = J {l + yi)dyi ...J (1 + y m )dy m ip N (xiy u . . .,x m y m ) e~ y i /2 ~"-~^ /2 (89) 



Another way to deal with divergence of resolvents, often used in practice, is to introduce the genus expansion: 
namely, to consider generating functions for fat graphs of fixed genus g: 



p {g) {x u ...,x m ) = £ KgL im x?...x% (90) 



i-i...i m —0 



where 



K it.i m = coefficient of N dcg(s) in K h ... im (N), deg(g) = (n + . . . + i m )/2 + (2 - 2g) - m (91) 



Recall, that genus g contribution to the connected correlator of tr c/f 1 . . . tr (f>' lm scales as N to the power 
(ii + . . . + «m)/2 + (2 — 2g) — m. Defined in this way, the genus g standard correlation functions (the genus 
g resolvents) are no longer divergent: the number of fat graphs of fixed genus grows much slower, than the 
total number of fat graphs. We will now demonstrate that relation (89) indeed allows to describe the genus 
expansion and find resolvents for any particular genus g. 



20 



6.2 1-point function 

The exact 1-point resolvent is given by 



oo 

Pn{x) = J dy e-y 2 / 2 <p N (xy) (92) 



where the 1-point Harer-Zagier function is given by 



2x 2 \\l-x 

Technically, to extract the genus expansion it is most convenient to introduce another variable X = x\/~N 
(see s. IV of [3]). As a consequence of the scaling rule (91), in terms of X the genus expansion becomes 
simply the l/N expansion. Indeed, in terms of X the Harer-Zagier function takes form 



X \ N (N + X 



y]v7 2x 2 l V n-x< 

and posesses an expansion in negative powers of TV: 

N + X 2 \ N f /N + X 2 \) / . 2A 6 2X W \ 2X 4k + 2 



N-X 2 ) ~ eXP Y'"*\N-X 2 J j 3N 2 ' 5N 4 ' ---J -^^( 2 jfc + i)JV2fc 

The exponent of a series can be represented as a series again: 

/ 2X6 2X io \ 2 / 2X 6 \x w + lx 12 lx u + ^-X^ + ^-X w 

exD 2X 2 + — + ^— + I - e 2X • 1 + — + 5 9 , 7 15 81 , 

cxp + 3N2 + ^ + ...J - e 11+ 37y2 + ^4 + ^6 + 

or, in a general form, 

/ ~ 2X 4k+2 \ _ ~P 1 X 4 P+ 2 g 2 g 

6XP Lrtl. , 1^21 " e 'l^l^ n \ M2p 1^ 



f , /N + X 2 \) / 2 2A 6 2A 10 \ 



Accordingly, the Harer-Zagier function takes form 



/A\ e 2X -1 e 2X X A e 2X fX 8 X w \ e 2X f X 12 ... 
<^ -== = A—-— + — - — + __ _ + _ +_- _ + __ + 



\y/N J 2X 2 N 3 A 3 \ 5 9 / A 6 \ 7 15 81 



or, generally, 



X \ _ e 2X ' 2 - 1 2X 2 ~ p 1 X^ 2 ^ 2 v 2"- 1 
VaJ ~ 2A2 +e "^Z^i A2p-i 2-. (2ii + 1) . . . (2i, + 1) ( } 

p=l 9=1 Ji + ...+t g =p x y ' 
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The leading contribution here is the genus Harer-Zagier correlation function: 



P(0)W = -2jp2- ( 95 ) 



Accordingly, the next-to- leading contribution corresponds to genus 1, and so on: 

= £ e "' (96) 

, . [X 12 2X 14 2X 16 \ 2X 2 , s 



Passing back to resolvents by taking Gaussian integrals (92), we obtain in genus zero 

OO o 9 

mm-J^^"^^- 1 -^^- <*») 

— OO 

This is the celebrated Wigner semi-circle distribution [20] . Similarly, in higher genera we obtain 

P(1) (X) =3X 4 (l-AX 2 y 5/2 (101) 
p (2 ){X) = (2lX 8 + 21X 10 )(l-4X 2 y 11/2 (102) 
p (3) {X) = (1485A 12 + 6138A 14 + 1738AT 16 ) (l - AX 2 ) ~ 17/2 (103) 



_ ^ (4g + 2g - 3)!!X 4 g+ 2 9- 2 (l - 4X 2 ) 1 / 2 - 2 g-" ^ 2^ 

P{9){ ) ~h *~ ' ~ " 



(104) 



As one can see, it is a straightforward exercise to extract the genus expansion from the Harer-Zagier 1-point 
function - it is even possible to write a formula for arbitrary g. Using the main result of present paper - the 
exact 2-point Harer-Zagier function - we can do a similar calculation at the 2-point level. 
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6.3 2-point function 

The exact 2-point odd and even resolvents are given by 



oo oo 

Pn{x,v) = J du J dv exp ^ - ^ iys+(xu,yu) (105) 



-oo — OO 



oo oo 



p N (x,y)= J udu J vdv exp ^ - ^ (p N (xu,yv) 



(106) 



-oo — oo 



where the 2-point odd and even Harer-Zagier functions are given by (56) and (57): 



¥>jv(&.y) = y 

o 



/ 

rfi / / 1 + t 2 + x 2 + y' 



2(x 2 +y 2 ) \\l + t 2 -x 2 -y 



1 



+ 7 syrff // l + t 2 + x 2 + j/ 2 \ W 2jV(x 2 +;/ 2 )(l+t 2 ) + (.x 2 + ? y 2 ) 2 -(l + t 2 ) 2 \ 

<Pn(*>V) J {x 2 +y 2 ) 2 \\ 1+t 2_ x 2_ y 2) (1+i 2 +a ,2 +y 2 )(1+ i2_. T 2_ y2) +M 

V 7 

The full resolvent is a sum of odd and even parts: 

Pn{x,v) = p { N ] {x,y) + p { x\x,y) (107) 



Again, to extract the genus expansion we introduce another variables A = x\/N, Y = yVN. As a conse- 
quence of the scaling rule (91), in terms of A and Y the genus expansion becomes simply the 1/N expansion. 
Taking the intermediate i-integral (we do not go into details here) we obtain in the odd case 

_ (_X_ J_\ _ XY(e 2X2+2Y2 -l) e 2X 2 +2Y 2 / X 3y3 x 5y X Y 5 \ 

LPn \Vn , Vn) ~ 2(X 2 + Y 2 ) + N 2 *v 3 + 3 + 3y + -- 



and in the even case 



2X 2 +2Y 2 / 4 4 4 OK 4 , n 4 



+ = — -x°y 2 + -A 4 y 4 + -x z y° + -A 4 r 2 + -A 2 y 4 + 

A 2 I 3 3 3 3 3 



We stop at genus one and do not write the genus two and higher contributions, it is absolutely straightforward 
to obtain them as higher 1/A corrections. To write a formula for arbitrary g (like we did in the 1-point case) 
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is a more involved, but still feasible exercise, which remains to be done. For the Harer-Zagier functions in 
genus zero, we obtain 



*U X > Y ) = { £Ty-y ( (2 * 2 + 2y2 - l)e2X2+2Y2 l ) (108) 



XY(e 2x2+2Y2 -l) 

^P(o)( x , Y ) = 2 (x 2 + y 2 ) 



(109) 



In genus one we have 



(110) 



Passing back to resolvents by taking Gaussian integrals (105) and (106), we obtain in genus zero 



+ X 2 F 2 / 1-2X 2 -2F 2 \ 

' j (x 2 -r 2 ) 2 {VT^X^VT^W^ J 



(112) 



P ( o)(*. y ) ~ 2(x2 _ y2)2 I Vl _ 4XV1 _ 4y2 - * y j ( 113 ) 



and in genus one 



+ (Y V\ - 8X2y4 + 8x4y2 + 8X 2 F 6 - 104X 4 F 4 + 8X 6 Y 2 - 32X A Y 6 - 32X 6 Y 4 + 6A0X 6 Y 6 
P(i)( X ' Y )- (1 - 4X 2 ) 7 /2(i _ 4y2)7/2 ( 114 ) 

- (l-4X 2 ) 7 /2(i_4y2)7/ 2 ( 115 ) 

Similarly, contribution of any higher genus can be found by expanding the exact functions (56) and (57). 
This method should be also applicable to 3-point and higher resolvents, but, unfortunately, 3-point and 
higher analogues of (56) and (57) are not found yet. 

7 Conclusion 

Despite the apparent simplicity and transparency of the Gaussian Hcrmitian model, which is beyond any 
doubt one of the most studied and best understood matrix models, its correlators form a complicated 
combinatorial system. Given a family of correlators, we can rarely explicitly describe its behaviour. Miriads 
of integer numbers, counting appropriate fat graphs or discrete Riemann surfaces, appear in a seemingly 
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random fashion. Integer numbers form patterns, they grow and they change according to laws which, despite 
the model is Gaussian, are far from being simple. In the case of one-point correlators 



tr cj) 2 jj = N 
tr =2N 3 + N 

tr cj) 6 ^ = 5N 4 + ION 2 
tr <f> 8 ^ = UN 5 + 70N 3 + 21N 
tr 6 10 )) = 427V 6 + 4207V 4 + 4837V 2 



these laws can be summarised in one compact formula, found by Harer and Zagier: 



tr (j) 2k > > A -j 



= coefficient of x 2k X N in 



(2fc-l)!! 1-A(l-A)-(1 + A):r 2 

The modest aim of our research was to find analogous formula for the two-point correlators 

((tr cj) tr ^ = N ((tr 4> tr 3 ^ = 37V 2 ^tr <f> tr </> 5 ^ = 107V 3 + 57V 

' tr <j)tr <j> 7 ^ = 357V 4 + 707V 2 
'tv , l.r <b 3 ^j = 127V 3 + 37V ^tr </) 3 tr = 457V 4 + 607V 2 

'tr 4? tr (f> 7 ^ = 1687V 5 + 6307V 3 + 1477V 
'tr </> 5 tr <j) 5 ) ) = 180N 5 + 6007V 3 + 1657V ^tr </> 5 tr </> 7 ^ = 7007V 6 + 49007V 4 + 47957V 2 
'tr d> 7 tr <b 7 )) = 28007V 7 + 343007V 5 + 813407V 3 + 166957V 



and it appears to be 



((tr^+ 1 tr^+ 1 \\ \ arctan ( — == |^=|= — 

11 = coefficient of x 2k+1 y 2m+1 \ N in ' VV ' ' 



(2i + l)!!(2j + l)!! " ' ' y (A-l) 3 / 2 - 1 + (A + l)(x 2 + y 2 ) 

This is of course just the first step (or, better to say, the second step). Three-point and higher correlators 
are still under-investigated. We are yet very far from complete understanding of integer numbers related to 
fatgraphs: hopefully, many more compact and beautiful formulas lie in wait. 
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